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Abstract 

We address the problem whether a given set of expectation values is 
compatible with the first and second moments of the generic spin operators 
of a system with total spin j. Those operators appear as the Stokes oper- 
ator in quantum optics, as well as the total angular momentum operators 
in the atomic ensemble literature. We link this problem to a particular 
extension problem for bipartite qubit states; this problem is closely re- 
lated to the symmetric extension problem that has recently drawn much 
attention in different contexts of the quantum information literature. We 
are able to provide operational, approximate solutions for every large spin 
numbers, and in fact the solution becomes exact in the limiting case of 
infinite spin numbers. Solutions for low spin numbers are formulated in 
terms of a hyperplane characterization, similar to entanglement witnesses, 
that can be efficiently solved with semidefinite programming. 

1 Introduction 

The "black magic" calculus of quantum mechanics allows us to make predictions 
about expectation values of certain measurement outcomes; however these ex- 
pectation values are surprising in the following ways: On one hand we know from 
Bell's inequality [1] that not all possible quantum mechanical expectation values 
are compatible with a local hidden variable model; on the other hand, not all 
expectation values which originate from a non-signaling constrained probability 
theory are quantum mechanical [2]. Given the quantum mechanical descrip- 
tion of the measurement device, the question arises which expectation values 
are compatible with quantum theory at all. The quantum mechanical moment 
problem [3, 4], as well as its truncated version [4] are famous paradigms of 
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this sort of question. An operational characterization of physical expectation 
values is of course desirable, however, only few examples are known so far. A 
well-known example is provided in the context of Gaussian states for systems 
with finite numbers of degrees of freedom. The first two moments of the cor- 
responding position and momentum operators are compatible with quantum 
mechanics if and only if the Schrodinger-Robertson uncertainty principle is ful- 
filled [5]. The condition can be written in terms of the covariance matrix [6], 
which allows — together with a vector formed by the mean values — an opera- 
tional low-dimensional description of all the quantum states compatible with 
the given moments. In most of the literature about Gaussian states this low- 
dimensional description is exploited, which pinpoints the importance of such a 
result. By contrast satisfying the Schrodinger-Robertson uncertainty principle 
is not sufficient to guarantee the existence of a quantum state in the case of 
higher moments [7], and a straightforward extension of Gaussian states might 
be non-trivial [8]. 

Clearly, given a set of linearly dependent operators the corresponding ex- 
pectation values have to reflect the same linear dependence. In addition, each 
expectation value has to be in the convex hull of the spectrum of the corre- 
sponding operator. In order to exploit further structures of the set of operators, 
we consider operator sets with an underlying Lie algebra structure. In the fol- 
lowing we consider operators which are irreducible representations of the Lie 
algebra su(2), the three-dimensional vector space of all tracefree, 2x2 complex 
matrices X with = —X. The corresponding irreducible representation on 
the Hilbert space Hj of dimension d = 2j + l is denoted by diTj : su(2) — > B(Hj). 
In addition, we use the label — diTj(o~k/2) for an irreducible representation 
of the Pauli operators o^, which constitute a basis for the Lie algebra su(2). 
This set of operators satisfies the commutation relations 



where sum denotes the Levi-Civita tensor. These commutation relations are 
well-known to be satisfied for the spin operators, and indeed the set of operators 
C^p can be considered as the familiar spin operators of a total spin j, and we 
use the term in the following. In fact, the considered scenario appears in a 
variety of different fields in physics: The macroscopic spin measurements of a 
state of an ensemble of N two-level atoms only supported on the symmetric 
subspace are spin operators of a spin j = N/2 system [9]; Any Stokes operator 
acting on a two mode system with fixed total photon number TV is described by 
a similar formalism [10], however in contrast to the spin operators the Stokes 
operators differ by a factor of 2 in the definition of the commutation relation, 
hence C[ N/2) =2C[ N/2) . Note, since there is only one irreducible representation 
of the Lie algebra su(2) in a given, fixed dimension d, the spin operators are 
unique up to unitary transformations. In addition, the spin operators satisfy 
the Casimir identity given by 
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In the following, wc arc interested in the expectation values of products of 
two spin operators only. Although higher moments can be measured in principle, 
in experiments it is often very tedious to get accurate values for those moments, 
e.g., Rcf. [11]. Let us formally denote the set of density operators p on a given 
Hilbert space H by 



The problem stated below asks for the compatibility of a given set of expectation 
values with spin operator measurements on a system of fixed dimension 1 . 

Problem 1.1 (Truncated su(2) moment problem). Consider the set of operators 
with k, I G {1, 2, 3} formed by the products of two spin operators acting 
on the Hilbert space Jij with dimension d = 2j + 1. Given a set of expectation 
values Mki S (D with k, I e {1, 2, 3}, under which conditions do these expectation 
values originate from a quantum mechanical state, i.e., 3 p € T>(Hj) such that 



We like to provide an operational description of the set of valid expectation 
values that enables working with the moments directly rather than using the 
complete density matrix. In order to check for compatibility of a given set of 
expectation values Mm, the following solution is provided: First, one verifies the 
linear dependence imposed by the Casimir identity. Next, one reconstructs a 
particular operator Pj{M) which acts on the symmetric subspace of two qubits. 
The given expectation values are consistent with the spin operators if and only 
if the operator Pj{M) represents a valid density operator for two qubits, that 
has exactly 2j — 2 Bosc-symmctric extensions. Although this reformulation does 
not provide an operational description yet, it opens the possibility to apply re- 
sults from this extension problem directly to Problem 1.1. Consequently, one 
can formulate operational approximations to the truncated moment problem. 
Whenever one finds a separable two qubit state pj(M), then it can be assured 
that the corresponding expectation values are quantum mechanical, while if one 
detects non-positivity of a certain operator Tj(pj(M)) ^ 0, that depends on the 
total spin number j and the reconstructed two-qubit state Pj(M), then the ex- 
pectation values are incompatible with quantum mechanics. This approximate 
characterization gets more accurate the larger the total spin number j becomes, 
and converges to the exact solution in the case of infinite spin numbers. 

The concept expectation value matrix is introduced in Sec. 2, which imposes 
already a strong condition on quantum mechanical expectation values in general. 
The Lie group structure simplifies the problem to a particular standard form 
of the given expectation value matrix. In Sec. 3, we relate Problem 1.1 to 
the characterization of Bose-symmetric extendible two qubit states by using 
a particular representation of the spin operators. In addition, we show how 
this idea directly provides an operational solution to a simplified version of the 

1 The considered set of given expectation values stands for knowing the first and second 
moments of the spin operators in arbitrary directions of the coordinate frame. 
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truncated moment problem. Since the exact solution to this extension problem 
is unknown, two different methods are considered in order to provide solutions 
for large (Sec. 4) and small spin numbers (Sec. 5). In particular, Sec. 4 deals 
with the two different approximation methods which both converge in the limit 
of high spin numbers. In Sec. 5, a solution to problems of the kind like Problem 
1.1 in terms of hyperplanes is provided, which can be efficiently solved for low 
spin numbers j. A graphical comparison between the different sets is given in 
Sec. 6 that demonstrates already the convergence. In Sec. 7 we summarize and 
give an outlook on possible further directions. 

2 Expectation value matrix 

Each expectation value Mki in Problem 1.1 is labeled by two indices, therefore 
each expectation value constitutes an entry of a particular matrix, which we term 
expectation value matrix in the following. Considering the set of expectation 
values M in matrix form already enables us to derive a strong statement about 
quantum mechanical expectation values in general; it even holds without the Lie 
group structure and has already appeared in the literature [12, 13, 14, 15, 16]. 
The construction of the expectation value matrix is equivalent to the derivation 
of the Schrodinger-Robertson uncertainty principle, which all given expectation 
values must clearly satisfy in order to originate from a quantum state. 

2.1 Definition & Properties 

The expectation value matrix defines a linear map on the set of density opera- 
tors, which preserves hermiticity and positive semidefiniteness. Therefore every 
expectation value matrix M, which originates from a valid quantum state p, i.e., 
M = M(p), must necessarily be positive semidefinite. The following proposi- 
tion introduces a slightly bigger expectation value matrix 2 \ '■ B(Hj) ~~ * B(<C 4 ), 
which contains our original expectation value matrix M in form of a submatrix, 
hence the given statements follows automatically. We skip the proof since it has 
already appeared in the literature, e.g., Ref. [13]. Note that further constraints 
are naturally imposed if the corresponding operator set shows linear depen- 
dence. Let us mention that Gaussian states are examples where the conditions 
on the expectation value matrix are necessary and sufficient for the existence of 
a quantum state. 

Proposition 2.1 (Expectation value matrix). Let T denote the set of operators 
acting on the Hilbert space Tij formed by the identity operator and the spin 
operators. The expectation value matrix \ '■ B{Hj) ~ * B(<C 4 ) is defined by 

B{Hj)3 A^ X ki(A)=tv(FlF l A) (4) 

with Fk G T for fc = 1, . . . , 4. This expectation value matrix \ has the following 
properties: 

2 We use the term expectation value matrix for both, the 3x3 matrix M and the 4x4 
matrix \. It should be clear form the context to which matrix one refers. 
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. X (aA + bB) = a X (A) + b X {B) for A, Be B(Hj), a,be<£. 

• IfA = A* then X (A) = X (AY . 

• If A > then \{A) > 0. 

2.2 Standard form 

In order to reduce the number of parameters in the expectation value matrix 
M we give a standard form first. The standard form allows us to consider ex- 
pectation value matrices with purely imaginary off-diagonal entries only, hence 
the set of expectation values is solely characterized by the mean values and 
the variances of the spin operators . In order to derive the standard form 
one has to show the covariance property of the expectation value matrix un- 
der the adjoint representation of the corresponding Lie group SU(2), given by 
Ad : SU(2) -» B(su(2)) and Adu(x) = WxU, for U G SU(2) and x e su(2). 
Acting on the (D 3 , this adjoint representation is unitary equivalent to the ir- 
reducible spin-1 representation m of the Lie group SU(2). This allows us to 
diagonalize the real part of the expectation value matrix. The idea is quite 
common in the literature of atomic ensembles, see, e.g., Ref. [17], hence the 
proof is omitted. 

Proposition 2.2 (Covariance property & Standard form). The expectation 
value matrix M — M(p) satisfies the covariance property 



where ttj denotes the irreducible representation of the group SU(2). Therefore 
each expectation value matrix M can be transformed to the standard form 



where D is a diagonal matrix with fixed trace, tv(D) = j(j + 1) and iA denotes 
the antihermitean, tracefree part. 

3 Reduction 

Each irreducible representation dnj of the Lie algebra su(2) can be built up 
from its two-dimensional, fundamental representations diri/z by acting on a re- 
stricted subspace of the tensor product space of multiple qubits [18]. Such a 
particular form of the representation appears as the usual representation in the 
case of atomic ensembles [9] . Given a valid density operator on this restricted 
tensor product space of multiple qubits, only an effective two qubit state suffices 
to determine the expectation values of the products of at most two Lie algebra 
elements, cf. Ref. [19]. Using the prescribed formalism enables us to establish a 
connection between Problem 1.1 and the Bose-symmetric extendibility question 
for a particular two qubit state, which constitutes the first result. This partic- 
ular extension problem shows similarities to the so-called symmetric extension 



M(nj{U)pnj(U)^) = 7ri(C/)M( / 9)7ri(J7) t , W e SU(2) 



(5) 



M = D + iA. 



(6) 
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problem, which recently has drawn attention in the literature [20, 21, 22, 23, 24]. 
Moreover, we further demonstrate the reduction idea by solving a simplified ver- 
sion of the truncated moment problem. 



3.1 Group theory 

This part reviews the irreducible representation dwj of the Lie algebra su(2) on 
the tensor product space of multiple qubits, restricted to the symmetric sub- 
space. The symmetric subspace H® 3 , indicated by the subscript +, consists 
of all vectors \ip) E H.® 2 ' J , with H = (D 2 , that arc invariant under all possible 
permutations TT2j(p)\ip) = \i ) ),P € &2j- Here tt2 j denotes an irreducible repre- 
sentation of the permutation group S2j , which acts on a basis of product states 
as n 2 j(p)\ii) ® ■ ■ ■ <g> \i2j) — Kp-i(i)) ® • • • ® \ip- 1 (2j))- The following proposition 
describes the particular form of the representation. 

Proposition 3.1 (Irreducible su(2) representation). One has: 

• The two Hilbert spaces (D 2j+1 and H® 2,3 with H = (D 2 are isomorph. 

• The irreducible representation of the Lie algebra dirj : su(2) — » B(H® 23 ) 
is given by 



with x( fc) = l® k 1 ® dni/2(x) ® l® 2 ^ fc , dni/2{x) is the two dimensional 
fundamental representation, and \ H ®2j denotes the restriction to H® 23 . 

3.2 Two qubit reduction 

Since the two abstract Hilbert spaces of the previous section are isomorph, any 
density operator of a spin j system can be interpreted as the density operator 
of a particular multipartite system of 2j qubits, i.e., pj e T>(Hf 23 ). In the 
following we are only interested in the expectation values of products of at most 
two spin operators. Because of the particular symmetry of the multipartite 
qubit state and the explicit representation of these operators, these expectation 
values depend only on an effective spin-1 system. The following proposition 
summarizes this reduction and links the corresponding spin operators. This 
reduction idea already appears in the context of atomic ensembles; a proof can 
be found, e.g., in Rcf. [19], which contains also examples of these reduced states. 

Proposition 3.2 (Spin-1 reduction). For every density operator lo G D(Tt® 23 ), 
the expectation values of the products of at most two Lie algebra elements only 
need to be calculated on an effective spin-1 state pj = tr3...2j(a>) € V(Hf 2 ) and 
are given by ti(diTj(x)u)) = tv(K.(d-Kj{x))pj) and similar tr(<9j7r(x)' 1 ' d-Kj{y)uj) = 




(7) 
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tr ■( A(d 'j it (x)^diTj(y))pj), for x,y e su(2) with the operators given by 3 



(8) 



for the mean values, and similar for the products of two spin operators 

A(dirj(x)tdnj(y)) = (2j dir 1/2 (x)idw 1/2 (y) ® 1 (9) 
+2j(2j - 1) a7r 1/2 (x)t ® a7T 1/2 (l/))) 

Only the reduced spin-1 state pj = tr 3 ... 2 j ;(w) of a given multipartite state u> 
of total spin j determines the expectation values of products of two spin oper- 
ators. However not all possible spin-1 density operators are reduced states of 
such multipartite states. The formal definition of valid two qubit reductions of 
a spin j system is given by 

Sj = { P e V(Hf )\3uj e V(Hf j ),tr 3 ... 2 » = p] , (10) 

with j > 1. Hence, a spin-1 state can correspond to a spin j system if and 
only if it can be extended to a multipartite state of 2j qubits which has sup- 
port only on the symmetric subspace. This definition resembles the notion of 
Bose- symmetric extensions for bipartite states, which is closely related to the 
symmetric extensions for bipartite states, cf. Ref [20] 4 . The next theorem 
establishes the connection between the Problem 1.1 and the problem of Bose- 
symmetric extensions of a particular spin-1 state. 

Theorem 3.1 (Bose-symmetric Extensions <=> Quantum mechanical expecta- 
tion values). Given the expectation value matrix M the corresponding spin-1 
state pj = Pj{M) is uniquely determined by 

tr(A(4 j) £p ) )p j )=M w , (11) 

for all k,l G {1,2,3} ; and A(-) is given by Eq. 9. The given expectation value 
matrix M is quantum mechanical if and only if Pj{M) e Sj, i.e., the state 
Pj{M) has exactly 2j — 2 Bose-symmetric extensions. 

3 Note, A(-) is used as a symbol here and should not be interpreted as a map A : £?(W® 2 "' ) — > 
B(H® 2 ). 

4 Note, that in the definition of symmetric extensions, as given in Ref. [20], the correspond- 
ing multipartite state must be invariant under the permutation of the individual subsystems, 
while in the case of Bose-symmetric extensions one requires for the multipartite state to have 
support on the symmetric subspace only. Although both problems are very similar — and in 
fact, many results can be "borrowed" directly form the symmetric extension case — certain 
properties differ. For example, every separable, permutation invariant density operator has a 
symmetric extension to arbitrary many copies, while it does not need to have a Bose-symmetric 
extension. 
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Proof. It is straightforward to check that the operators given by Eq. 9 span an 
operator basis for a generic spin-1 state, hence the corresponding spin-1 state 
Pj = Pj{M) is completely determined by the expectation value matrix M. This 
particular two-qubit state can correspond to a spin j system if and only if it is 
an element of the class Sj . □ 



3.3 First moment problem 

In the present subsection, we briefly discuss a simplified version of the truncated 
moment problem in which one asks for compatibility of the first moments only. 
This problem allows an operational solution and highlights again the importance 
of the reduction idea. It states as follows: 

Problem 3.1 (First moment problem). Consider the set of spin operators L k ^ 
with k € {1, 2, 3} acting on the Hilbert space Tij of dimension d = 2j + 1. Given 
a set of expectation values Lk € 1R with k G {1,2,3} ; under which conditions 
do these expectation values originate from a quantum mechanical state, i.e., 
3pe V{Hj) such that tr{C k j) p) = L k ,yk? 

The operational description of the set of valid expectation values relies on 
the reduction to a spin-1/2 problem, in analogy to Proposition 3.2, and that 
every valid spin-1/2 state represents a possible reduced state of a particular 
spin-j state, i.e., each single qubit density operator can be extended to a Bosc- 
symmetric density operator of 2j qubits. The next proposition contains the 
operational solution. 

Proposition 3.3 (Operational description). The expectation values L k € R 
with k e {1,2,3} are compatible with the spin operators of a generic spin j 
system if and only if it holds 

Ei?<J 2 (12) 

Proof. In analogy to Proposition 3.2, it holds that for every spin-j density op- 
erator u> e V(Ti.f 2: '), the expectation values of each Lie algebra element d-Kj{x) 
need only to be calculated on an effective qubit state pj = tr2...2j(w) and are 
given by tv(dTTj(x)uj) = tv{k{d-Kj{x))p~j) with A(x) = 2jd r K\/2{x). Contrary to 
the spin-1 case, every qubit state p can be extended back to a multipartite state 
u) = p® °\-,®2j, i.e., to a density operator of a spin j system with the correct 

first moments. Thus, the given expectation values L k only have to form a valid 
qubit density operator Pj({Lk}) € V(TL). The set of given expectation values 
Lk uniquely determines the Bloch vector Lk/ j of the reduced qubit state pj, 
which is quantum mechanical if and only if Eq. 12 holds. □ 

We close with a discussion about possible solution states to the first moment 
problem. Given the set of expectation values, any compatible quantum state can 
be parameterized as p(x) = ps^{{Lk}) + /O op on(x), with a fixed part pfi x ({£fe}) = 
l/(2j + 1) + ^2 k o-kC^k an d = Lk/\\£h tnat i s completely determined by 
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the first moments Lk, and an orthogonal open part p pen(x), which linearly 
depends on some open parameters x, Xi € R, and which have to be chosen 
such that p(x) > forms a valid density operator. Proposition 3.3 guarantees 
existence of such a set of open parameters x as long as the given expectation 
values fulfill Eq. 12. A special class of solutions constitutes the case in which 
PRx({Lk}) > forms already a positive semidefinite operator itself. However, 
this is the case if and only if 5Z i=1 Lf < (j + l) 2 /9, which shows that the open 
part Popon(x) is indeed necessary since it allows compensation of a non-positive 
fixed part pfi x ({ifc}) ^ in certain cases. 



4 Approximating sub- and superset 

Any solution to the Bose-symmetric extension problem for two qubits consti- 
tutes a solution for the su(2) moment problem; however, an analytic character- 
ization of those sets appears cumbersome. Therefore, two operational approxi- 
mation methods to the set of extendible states are described in the following, 
which allows identification of certain expectation value matrices, that are either 
quantum mechanical (subset) or not (superset). Both approximating sets con- 
verge to the exact set in the limit of infinite numbers of extensions. The subset 
characterization relies on a result by Doherty et at, Ref. [20], while the super- 
sets depend on the expectation value matrix. Both approximation sets (and 
the corresponding convergence statements) are given in terms of renormalized 
expectation values of an expectation value matrix defined by 

_</#>) _ 1 

Uk = — Vk = liTW)~^-i' (13) 

for k G {1,2,3}. Note that one has only to consider these 6 expectation values 
because of the standard form of the expectation value matrix as described by 
Eq. 6. More precisely, the parameters Vk determine the diagonal entries of 
the expectation value, while Uk fixes the antihermitean, tracefree part. These 
renormalized expectation values have the advantage that the determined spin-1 
state pj({uk, Vk}) is independent of the considered total spin number j. 



4.1 Approximating subset 

The set of separable two qubit states, supported on the symmetric subspace, 
forms an approximating subset to the set of Bose-symmetric extendible states, 
independent of the number of extensions. Hence, if a given expectation value 
matrix M allows the reconstruction of a separable spin-1 state pj(M), then it 
can be assured that all expectation values in M are quantum mechanical. The 
set of separable spin-1 states is given by 

n = [ P G V{Hf)\ P T > o} , (14) 

where T denotes the partial transpose, which is necessary and sufficient to char- 
acterize separable states in these low dimensions [25], [26]. The next theorem 
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outlines the properties of this subset. The convergence result is a direct conse- 
quence of a corresponding result for states which have a symmetric extension 
[20]. For a fixed number of photons, any convex combination of polarization 
states, where all photons are in one optical mode, constitutes physical relevant 
examples where the qubit density operator is separable. In contrast, any spin 
squeezed state must have an entangled two qubit density operator [27], as well 
as certain multiqubit states in the atomic ensemble literature [28]. 

Theorem 4.1 (Approximating subset). The sets Sj with j > 1, and TZ given 
by Eqs. 10 and 14 respectively satisfy: 

1. If p G TZ, then p G Sj,Vj. Hence TZ C <Sj,Vj. 

2. If p G Sj, then p G Sj' with j' < j. Hence Sj C Sj' for j' < j. 

3. lim^oo^ = f]Sj C TZ. 

j 

Proof. Any separable density operator p G TZ can be written as a convex 
combination of pure product states. We assume the following decomposition 
p = ^2iPi\cti){ai\ ® \Pi)(Pi\, where {pi} denotes the probability distribution and 
\a.i), G H,Vi. Because of p e B{Hf 2 ), it must hold that tr(P_p) = 0, where 
P_ > denotes the projector onto the antisymmetric subspace. Given the 
particular decomposition into positive semidefinite operators it follows that the 
trace must vanish for each term, tr(P_ £g> = 0,Vi, which shows 

that each |ai)(g)|/3i) G H% 2 . This can be the case iff \cti) = |/?i),Vi. This allows us 
to write down a valid Bose-symmetric extension given by oj = ^Pi|aj)(a;j|® 2: ', 
which then proves the first statement, p G Sj. 

The second statement follows trivially. The last point is a direct consequence 
of Theorem 1 in Ref. [20] , which states that any bipartite mixed states pab G 
T>(TLa ® Hb), which has arbitrary many symmetric extensions to one of the 
subsystems must necessarily be separable. □ 

4.2 Approximating superset 

Exploiting the characteristics of the expectation value matrix allows us to for- 
mulate an approximating superset to the set of extendible states. Any reduced 
spin-1 state pj = tr 3 ... 2 j of a given state uj G V{H.^) allows the computa- 
tion of the expectation values of a corresponding spin j system via the operator 
sets A(-) given by Eqs. 8, 9; in this sense the expectation values from a spin j 
system can be "recovered" from the reduced two qubit state. Proposition 4.1 
introduces the reduced expectation value matrix Tj, which tries to reproduce a 
normalized version of the original expectation value matrix \ of the spin j state 
only from the reduced spin-1 state; renormalization protects against divergence 
of the expectation values in the limit j ^ oo. It holds that if the two qubit 
density operator p has a valid extension to 2j qubits, then one obtains a positive 
semidefinite operator. 
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Proposition 4.1 (Reduced expectation value matrix). Let T denote the set of 
operators acting on the Hilbert space rij formed by the identity and the renor- 
malized spin operators / j . The reduced expectation value matrix of order 
j, denoted by Tj : B(Tif 2 ) — > £>((D 4 ), is defined as 



with Fk € T for k — 1, ... 4, and the operators A(-) are given by Eqs. 8, 9. The 
map is linear and preserves hermiticity. In addition, it holds that if p € Sj then 



The approximating superset consists of all spin-1 states for which the reduced 
expectation value matrix of order j delivers a positive semidefinite operator, and 
is denoted by 



for j > 1. Unlike the subset characterization, any superset depends on the 
number of considered extensions j. The next theorem lists the properties of 
the approximating superset. If a given reconstructed spin-1 operator Pj(M) is 
not an element of the class Tj, then the corresponding expectation values in 
M are incompatible with quantum mechanics. The supersets converge to the 
set of separable states in the limit of infinite number of extensions, hence the 
reduced expectation value matrix delivers an alternative necessary and sufficient 
entanglement criterion, valid only for symmetric two qubit state, but which does 
not rely on partial transposition. However, in order to prove convergence, we 
exploit a necessary and sufficient criterion for entanglement of the reduced two 
qubit state in terms of the expectation values of the spin operators, cf. Ref. [29], 
which is derived using the partial transposition. 



Theorem 4.2 (Approximating superset). The sets Tj with j > 1 satisfy: 

1. If p <G 7), then p e Ty with j' < j. Hence Tj C Tj> for j' < j. 

2. If pe Sj, then p G Tj,\/j. Hence Sj C Tj,\/j. 

3. lim^oo T 3 ■ = f| T 3 ■ C K. 



Proof. The proof of the first two statements is straightforward. Convergence 
is proven via contradiction along the following line: First, one uses the result 
given in Ref. [29] to provide a necessary and sufficient criterion for two qubit 
entanglement in terms of the renormalized expectation values. Second, one 
derives a necessary condition for non-negativity of the reduced expectation value 
map Tj in the limit j — > oo, expressed again in terms of the renormalized 
expectation values of the spin-1 state. The two given conditions are mutually 
exclusive, hence no entangled state can be part of the superset Tj in the limit 
j — ► oo. This shows convergence to the separable states in the end, since 
separable spin-1 states are trivially part of the outer approximation. 



B(Hf 2 ) 3iH Tj (A) kl = tr(A(^ Fi)j4 ), 



(15) 



T» > 0. 



(16) 




(17) 
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For a given direction n € R 3 the corresponding spin operator in this direction 
is defined by = d r Kj{h ■ a/2), and a denotes the vector of Pauli matrices. 
Let Wft and v„ denote the corresponding renormalized expectation values of this 
spin operator given by Eq. 13. Using the result given in Rcf. [29], the reduced 
two-qubit state Pj only supported on the symmetric subspace is entangled iff 
there exists a direction h, such that — u\ < holds. 

The reduced expectation value matrix of order j of a given spin-1 density 
operator p, Tj(p) is positive semidefmite iff all principle minors are non- negative. 
Using the unitary freedom of the matrix representation, in combination with the 
non-negativity of all possible 2x2 submatrices ensures that Var(£l i) )/.?' 2 > 
for all possible directions n. If one re-expresses this condition in the normalized 
expectation values one arrives at 

«fl-«ft + ^:(l-«ft)>0, (18) 

which becomes — u\ > in the limit j — > oo. This condition ensures that 
the entanglement condition can never be met in the limit, and this proves the 
convergence. □ 



5 Hyperplane characterization 

In this part we discuss some general features of compatibility problems, in which 
one asks whether a certain set of expectation values are consistent with an 
operator set, similar to Problem 1.1. Consider a set of linear operators on a 
finite dimensional Hilbert space. Without loss of generality, only a set of linear 
independent, hermitian operators needs to be considered, because any operator 
can be decomposed into hermitian operators, and linear dependencies in the 
operator set demand only trivial conditions on the corresponding expectation 
values. 

Problem 5.1 (Consistency problem). Consider a set of linearly independent, 
hermitian operators acting on a finite dimensional Hilbert space H, {Ai} with 
i = 1, . . . , n. Given an expectation value vector b € R™, under which conditions 
do these expectation values originate from a quantum mechanical state, i.e., 
3pe V(H) such that tr{A iP ) = bi,Vi? 

Since the set of density operators V(TL) is compact, and because of linearity 
of the map Ai : T>(H) — > R™ that assigns the corresponding expectation values 
to a given state, p M(p)i = tr(Aip), the set of valid expectation values 
is also compact. According to the Hahn-Banach theorem, every expectation 
value vector outside this set can be separated by a corresponding hyperplane, 
in analogy to entanglement witnesses [26] for the separability problem. The 
following theorem summarizes the necessary properties of those hyperplanes and 
therefore provides a possible characterization of the set of quantum mechanical 
expectation values. Given a particular expectation value vector, the search for 
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the optimal hypcrplane can be cast in the form of a semidefinite program [30] , 
which can be solved efficiently with interior-point methods. The given proof is 
based on this idea and employs well-known duality relations from semidefinite 
programming, cf. Ref. [30]. 

In the case of the considered su(2) moment problem, a formulation in terms 
of a semidefinite program offers an efficient way to characterize the set of possi- 
ble expectation values for low total spin numbers j, i.e., in situations where the 
approximation method fails to provide a good description. Note that if one uses 
the particular representation of the considered spin operators, then one recov- 
ers the characterization of Bose-extendible states in terms of witness operators 
similar as described in Ref. [20]. 

Theorem 5.1 (Hyperplane characterization). Given a set of linearly indepen- 
dent operators {Ai} and an expectation value vector b <E R™, it holds: 

• The vector b is non-quantum mechanical, if and only if there exists a 
hyperplane, characterized by the normal vector z e R n : Z = Y] ZiAi > 0, 
tr(Z) = 1, which detects it z T h < 0. 

• The vector b is quantum mechanical, if and only if for all hyperplanes, 
characterized by the normal vector z <G R™ ■ Z z iAi > 0, tr(Z) = 1 it 
holds that z T h > 0. 

Proof. The Gram-Schmidt orthogonalization allows to transform any set of arbi- 
trary linear independent, hermitian operators {Ai} to an orthonormal set of her- 
mitian operators {S^}, i.e., the operators satisfy tr(SiSj) — Sij and tr(Si) oc Sn. 
The connection is 1-to-l and the corresponding expectation value vector b has 
to be transformed in the same way, which results in the new expectation value 
vector t. 

The set {Si} for i = l,...n can be extended to a hermitian operator ba- 
sis {Si} with i = l,...d 2 and d = dim(7Y). This basis set {Si} satisfies 
again tr(SiSj) = Sij and ti(Si) cx Sn, which states that the first operator 
is proportional to the identity, Si = 1/Vd, and all other elements Sj for 
j > 1 are tracefree. Therefore every density operator can be expressed as 
p = J2i x iSi with Xi — ti(Sip). It must hold Xi = U for i = 1, . . . n since oth- 
erwise, the expectation values do not match. The remaining open parameters 
Xi , Vi = n + 1 , . . . d 2 must be chosen such that p(x) = pa x + p opcn (x) > forms a 
positive semidefinite operator and we have defined the fixed part pfi x = J^^iSi 
and open part p pcn(x) = ^2xiSi of the density operator. This can be cast 
into the form of a semidefinite program [30], given by minuet subjected to 
F(x, t) = p{x) +tt > 0, with solution (i*,x*). If t* > 0, then there exists no 
parameters x such that p(x) > 0, since otherwise t* is not optimal. On the 
other hand, if t* < then p(x*) > 0, and the expectation value vector t is 
quantum mechanical. 

Every semidefinite program has an associated dual program, cf. Ref. [30], 
which reads maxz(- tr(Zpfi x )) subjected to Z — Y^ii=i z i^i > and tr(Z) = 1. 
Using the orthogonality of the operators Si, the objective value can be writ- 
ten as tr(Zpfi x ) = z T t, and its optimal value is denoted by d* . Note that both 



13 



semidefinite programs are strictly feasible: If one selects t > | min A(pfi x )|, where 
A(-) denotes the corresponding eigenvalues, one obtains a strictly positive solu- 
tion F(x = 0, t) > for the first program; the operator Z = t/d > provides a 
strictly positive solution for the second, dual program. These conditions ensure 
strong duality, which states equality of both programs t* = d*, and comple- 
mentary slackness, which guarantees that there exist actual parameter sets for 
(x op t,i pt), and Z opt that attain the solutions t* and d* respectively, cf. Ref. 
[30] . Therefore if t is not quantum mechanical it holds that 

-t* = -d* = tr(Z opt p fix ) = z£ t t < 0, (19) 

Hence every non-quantum mechanical expectation value vector is detected by 
the corresponding hypcrplanc. In contrast, if t is quantum mechanical, the 
weak duality condition [30] ensures that for every feasible solution z of the 
second program one has 

z T t = tr(Zp fix ) > -t* > 0, (20) 

so no quantum mechanical expectation value vector is detected by the hyper- 
plane. Using the described 1-to-l correspondence allows to translate these con- 
ditions back to the original operator set {Ai}, which delivers the result. □ 

6 Visualization 

The sets of renormalized expectation values that are determined by the de- 
scribed inner and outer approximation enclose the set of quantum mechanical 
expectation values, and convergence is reached in the case of infinite spin num- 
bers j. In order to visualize each convex sets, one starts with given expectation 
values Uk, k e {1,2,3}, which fulfill the condition ^2 k u 1 — 1> which ensures 
existence of a quantum state, cf. Proposition 3.3. Next, the remaining open 
parameters Vk must be chosen according to the conditions from either an ap- 
proximation set or from the exact solution. Using the Casimir identity, given by 
Eq. 2, allows to further reduce the number of considered parameters Vk, and we 
choose V\ and V2 in the following. For both approximations, the search for the 
extremal values of the parameters can be cast into the form of a semidefinite 
program, which can be efficiently solved by standard semidefinite program mod- 
ules [31], [32]; for the exact solutions one needs to employ a semidefinite program 
anyway, cf. Sec. 5. Figure 1 shows the exact and the approximation sets for 
the case of u = (0.1,0.2,0.3) and total spin number j = 5. The inset figures 
visualize the same sets for different spin numbers j, and should demonstrate 
the corresponding convergence. The outer approximation seems to describe the 
actual set of quantum mechanical expectation values with increasing accuracy 
as the total spin number j becomes larger. 
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Figure 1: Sets of renormalized expectation values v\ and v 2 described by 
the approximation in comparison with the exact solution for the case of 
u = (0.1,0.2,0.3). The inset figures are drawn with the same axis scale, and 
should demonstrate the effect of different total spin numbers j on the different 
sets. 



7 Conclusion 

We have addressed the problem whether a given set of expectation values can be 
compatible with expectation values for products of two spin operators. Those 
operators, as abstractly introduced as the irreducible representations of an un- 
derlying su(2) Lie algebra, appear for example as the Stokes operators in the 
quantum optics literature or as the total angular momentum operators in the 
case of atomic ensembles. Because of the particular product structure one can 
already impose a strong conditions on the given set of expectation values, which 
is summarized in positive semidefinitcness of a corresponding expectation value 
matrix that ensures the Schrodingcr-Robcrtson uncertainty principle. Exploit- 
ing a particular representation of the su(2) Lie algebra, allows us to relate the 
Problem 1.1 to the Bose-symmetric extension problem for qubits, hence the 
following solution to Problem 1.1 is provided: Suppose that a given set of ex- 
pectation values M satisfies the linear equality imposed by the Casimir identity, 
one reconstructs a particular two qubit operator pj(M). This operator pj(M) 
must represent a valid two qubit density operator of a bipartite qubit system 
with at least 2j — 2 Bose-symmetric extensions, otherwise the given expecta- 
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tion values M disagree with the predictions of quantum mechanics. Since an 
exact characterization of two qubit states with a definite number of extensions 
is cumbersome, we consider in particular the two extreme cases. In particu- 
lar, for large spin numbers we have presented two different approximating sets. 
Whenever one finds a separable two qubit state pj(M), then the expectation 
values can be assured to be quantum mechanical. Contrary, if one finds non- 
positivity of a particular expectation value matrix Tj(pj(M)) ^ 0, that depends 
on the total spin number j and the reconstructed two-qubit state pj(M), then 
the corresponding expectation values are incompatible with the spin operators. 
In combination, both tools allow an approximate operational description of the 
quantum mechanical expectation values. In particular, the presented method 
gets better the larger the total spin number j becomes, and convergence is as- 
sured in the limit of infinite numbers of extensions. In order to provide a feasible 
solution for low spin numbers, we characterize the sets of physical expectation 
values, similar as in Problem 1.1, via hyperplanes. The search for the optimal 
hyperplane can be cast into the form of a semidcfinite program which can be 
solved efficiently. 

It remains open whether one can find an operational characterization of Bose- 
symmetric extendible two qubit states. In addition, it might be interesting to 
further investigate whether similar ideas can be used if one considers different 
Lie groups; we leave this open for future discussions. 
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